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ABSTRACT 

The  Hall  thruster  cluster  is  an  attractive  propulsion  approach  for  spacecraft 
requiring  very  high-power  electric  propulsion  systems.  This  article  presents  plasma 
density,  electron  temperature,  and  plasma  potential  data  collected  with  a  combination  of 
triple  Langmuir  probes  and  floating  emissive  probes  in  the  plume  of  a  low-power,  four- 
engine  Hall  thruster  cluster.  Simple  analytical  formulas  are  introduced  that  allow  these 
quantities  to  be  predicted  downstream  of  a  cluster  based  solely  on  the  known  plume 
properties  of  a  single  thruster.  Ion  energy  distribution  functions  measured  using  both  a 
parallel  plate  electrostatic  analyzer  and  a  retarding  potential  analyzer  are  presented.  A 
cluster  of  Hall  thrusters  is  shown  to  exhibit  dramatically  different  ion  energy  profiles 
compared  to  a  single  thruster.  In  particular,  clustering  causes  a  significant  increase  in  the 
fraction  of  ions  at  energies  below  the  primary  peak  in  the  distribution,  most  likely  due  to 
an  increase  in  elastic  scattering  and  the  effects  of  the  unique  plasma  potential  profiles  in 
the  cluster  plume. 


Approved  for  public  release;  distribution  unlimited. 
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Nomenclature 


A  =  Area  of  one  electrode 

B  =  Magnetic  field  strength 

d  =  Energy  analyzer  plate  separation  distance 

e  =  Electron  charge 

E  =  Electric  field  strength 

kb  =  Boltzmann’s  constant 

K45  =  Spectrometer  constant  for  45°  energy  analyzer 

L  =  Distance  between  slits  in  energy  analyzer  baseplate 

mi  =  Ion  mass 

me  =  Electron  mass 

ne  =  Electron  number  density 

no  =  Reference  density 

qi  =  Ion  charge 

Te  =  Electron  temperature 

Ui  =  Initial  ion  velocity 

Vb=  Ion  beam  voltage 

Vd2  =  Voltage  measured  between  triple  probe  electrodes  1  and  2 
Va  =  Voltage  applied  between  triple  probe  electrodes  1  and  3 
V  =  Equivalent  ion  voltage,  miUi2/(2qi) 

Vf  =  Floating  potential 

VR  =  Energy  analyzer  reflector  plate  voltage 


<j)  =  Plasma  potential 
4>t  =  Thermalized  potential 

j  —  Subscript  denoting  the  contribution  from  an  individual  thruster 

Introduction 

Future  space  missions  will  require  electric  propulsion  systems  capable  of 
operating  at  very  high  power  levels  compared  to  those  currently  in  use. 1,2  One  method 
being  considered  for  reaching  these  power  levels  involves  clustering  multiple 
moderately-powered  devices  together  to  reach  the  total  throughput  desired.  The  most 
viable  type  of  electric  propulsion  device  for  this  class  of  mission  is  the  Hall  thruster  due 
to  its  low  specific  mass,  high  thrust  density,  and  high  reliability.1,2  In  an  effort  to 
understand  the  technical  issues  related  to  operating  multiple  Hall  thrusters  in  close 

proximity  to  each  other,  a  cluster  of  four  Busek  BHT-200-X3  200- watt  class  devices  is 
being  studied.3"5 

A  cluster  of  thrusters  may  have  a  slightly  lower  efficiency  and  higher  dry  mass 
than  a  single,  similarly  powered  thruster  since  larger  engines  have  historically 
outperformed  smaller  thrusters.  A  cluster,  however,  has  several  advantages  over  a 
monolithic  thruster,  including  improved  system  reliability  due  to  the  inherent  redundancy 
of  running  multiple  engines  and  the  ability  to  throttle  the  system  by  simply  turning  off 
one  or  more  thrusters.  Throttling  the  system  in  this  way  allows  the  cluster  to  operate  at 
lower  power  without  running  any  of  the  individual  thrusters  at  off-design  conditions. 

This  characteristic  of  a  cluster  may  prove  beneficial  on  missions  where  either  the 
available  power  or  the  propulsive  needs  change  as  a  function  of  time.  For  example,  a 
high-power  cluster  of  Hall  thrusters  could  be  used  for  the  initial  low-earth  orbit  to 
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geosynchronous  orbit  (LEO-GEO)  transfer  of  a  geosynchronous  communications 
satellite.  Upon  reaching  its  final  destination,  one  element  of  the  cluster  could  then  be 
used  for  north-south  station  keeping.  A  final  advantage  of  clustering  is  the  high  degree  of 
system  scalability.  In  principle,  once  the  technical  issues  involved  with  operating  a 
cluster  are  fully  understood,  a  single  flight-qualified  engine  could  support  a  wide  range  of 
missions  requiring  various  power  levels  by  simply  clustering  the  appropriate  number  of 
thrusters.  Thus,  enhanced  scalability  and  flexibility  make  clusters  attractive  for  many 
missions. 

Although  using  a  cluster  of  commercially  available  thrusters  for  primary 
propulsion  appears  to  be  advantageous  for  some  missions,  there  are  several  systems 
integration  issues  that  must  be  addressed  before  clusters  can  be  used  in  flight.1,2  In 
particular,  it  is  imperative  that  the  interaction  of  the  plasma  plumes  both  among  the 
thrusters  and  with  the  spacecraft  be  understood.  In  an  effort  to  address  this  issue,  the  ion 
energy  distribution  downstream  of  a  low-power  Hall  thruster  cluster  is  studied  using  both 
a  parallel-plate  electrostatic  analyzer  (ESA)  and  a  retarding  potential  analyzer  (RPA). 
Additionally,  the  electron  number  density,  electron  temperature,  and  plasma  potential 
downstream  of  the  cluster  are  measured  using  a  combination  of  electrostatic  probes.  In 
each  case,  the  profiles  recorded  in  the  cluster  plume  are  compared  to  those  measured 
downstream  of  an  individual  thruster. 

Experimental  Apparatus 

Cluster 
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The  cluster  used  in  this  experiment  is  composed  of  four  Busek  BHT-200-X3  200- 
watt  class  Hall  thrusters.  An  earlier  version  of  this  thruster  is  reported  to  operate  at  an 
anode  efficiency  of  42%  and  specific  impulse  of  1300  seconds  while  providing  12.4  mN 
of  thrust  at  the  nominal  operating  conditions.6  Each  thruster  has  a  mean  diameter  of  21 
mm  and  is  operated  on  xenon  propellant.  The  thrusters  are  arranged  in  a  2x2  grid  with 
approximately  1 1.4  centimeters  between  the  centerlines  of  adjacent  thrusters.  Typical 
operating  conditions  for  the  BHT-200  are  given  in  Table  l.3  Figure  1  shows  the  cluster 
during  operation. 


Table  1 :  Typical  operating  conditions  for  the  BHT-200  Hall  thruster. 


Parameter 

Value 

Discharge  Voltage  (V) 

250  ±0.5 

Discharge  Current  (A) 

0.80  ±  0.03 

Cathode  Potential  (V) 

-8.5  ±1.0 

Electromagnet  Current  (A) 

1.0  ±0.03 

Keeper  Current  (A) 

0,5  ±  0.05 

Keeper  Voltage  (V) 

13  ±  1 

Anode  Mass  Flow  Rate  (seem) 

8.5  ±0.85 

Cathode  Mass  Flow  Rate  (seem) 

1.0  ±0.1 
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Vacuum  Facilities 


Fig,  1 :  A  low-power  cluster  in  operation. 


The  electron  number  density,  electron  temperature,  and  plasma  potential 
measurements  presented  in  this  paper  are  obtained  in  Chamber  6  at  the  Air  Force 
Research  Laboratory  (AFRL).  Chamber  6  is  a  1.5  x  2.4  meter  cylindrical,  stainless  steel 
vacuum  chamber  that  is  evacuated  by  one  dual-stage  cryopump  and  four  single-stage 
cryopanels.  During  thruster  operation,  the  chamber  pressure  stabilizes  at  approximately 
6. lxl O'6  Torr  for  single  thruster  operation  and  2.3xl0‘5  Torr  for  four-thruster  operation. 
Both  reported  pressures  are  corrected  for  xenon. 

Measurements  of  ion  energy  distributions  are  conducted  in  the  Large  Vacuum 
Test  Facility  (LVTF)  at  the  University  of  Michigan’s  Plasmadynamics  and  Electric 
Propulsion  Laboratory  (PEPL).  The  LVTF  is  a  6x9  meter,  cylindrical,  stainless  steel  clad 
vacuum  chamber  that  is  evacuated  by  seven  cryopumps,  which  provide  a  pumping  speed 
of  500,000  liters  per  second  on  air  and  240,000  liters  per  second  on  xenon  for  typical 
base  pressures  of  approximately  2.5x10 7  Torr.  For  these  experiments,  only  four 
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cryopumps  were  used  resulting  in  chamber  pressures  of  1 , lxl O'6  and  3. 6x1  O'6  Torr 
(corrected  for  xenon)  during  single-  and  four-thruster  operation,  respectively. 

Coordinate  System 

The  naming  convention  and  coordinate  system  used  throughout  this  experiment 
are  shown  in  Fig.  2.  As  shown,  the  thrusters  are  labeled  as  TH  1-4  beginning  in  the  upper 
left-hand  comer  and  proceeding  counterclockwise.  The  origin  of  the  coordinate  system  is 
defined  as  the  midpoint  of  the  cluster  in  the  displayed  X-Y  plane.  The  Z  coordinate 
measures  the  distance  downstream  of  the  thruster  exit  plane.  A  three-dimensional 
positioning  system  is  used  to  sweep  probes  through  the  plasma  plume. 


Fig.  2.  Thruster  naming  convention  and  coordinate  system. 

Triple  Probe 

The  triple  Langmuir  probe  used  for  these  experiments  consists  of  3  tungsten 
electrodes  insulated  from  each  other  by  an  alumina  rod.  Each  electrode  is  0.5  mm 
(0.020  )  in  diameter  and  5.0  mm  (0.20”)  long.  The  spacing  between  the  centerlines  of 
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adjacent  electrodes  is  approximately  2  nun.  The  probe  is  sized  to  criteria  that  allow  the 
standard  “thin  sheath”  assumptions  of  probe  theory  to  be  applied.7  These  criteria,  which 
are  discussed  elsewhere,8  are  necessary  to  ensure  proper  operation  of  the  probe. 

The  symmetric  triple  probe,  originally  developed  by  Chen  and  Sekiguchi,9  is  a 
convenient  plasma  diagnostic  for  collecting  large  amounts  of  data  due  to  the  elimination 
of  the  voltage  sweep  required  by  other  electrostatic  probes.  Additionally,  since  the  probe 
as  a  whole  floats,  the  disturbance  to  the  ambient  plasma  is  minimized  compared  to  single 
Langmuir  probes,  which  draw  a  net  current  from  the  discharge.  The  relations  used  to 
determine  plasma  parameters  from  the  measured  quantities  are  given  in  Eqns.  1  and  2. 
Various  error  analyses  indicate  that  the  uncertainty  in  the  calculated  electron  temperature 
and  number  density  are  generally  less  than  30%  and  60%,  respectively.9, 10  The  relative 
uncertainty  between  multiple  data  points  recorded  using  the  same  probe  is  believed  to  be 
significantly  lower  than  the  absolute  uncertainty. 
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Emissive  Probe 

Plasma  potential  measurements  are  conducted  using  a  floating  emissive  probe 
similar  to  the  one  described  by  Haas  et  a/.11  The  emitting  portion  of  the  probe  consists  of 


a  0.127  mm  (0.005”)  diameter  tungsten  filament  loop,  the  ends  of  which  are  inserted  into 
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double  bore  alumina  tubing  along  with  0.508  mm  (0.020”)  diameter  molybdenum  wire 
leads.  Short  lengths  of  tungsten  wire  are  inserted  into  the  alumina  tube  to  insure  contact 
between  the  emitting  filament  and  molybdenum  leads.  The  diameter  of  the  emitting 
filament  loop  is  approximately  3  mm.  Figure  3  shows  a  sketch  of  the  emissive  probe. 

Molybdenum  Leads 


Alumina  Insulator 

Fig.  3:  The  floating  emissive  probe. 

The  theory  of  the  emissive  probe  is  well  established  and  results  in  the  conclusion 
that  a  thermionically  emitting  filament  will  assume  the  local  plasma  potential  when  its 
emitted  electron  current  is  sufficient  to  neutralize  the  plasma  sheath.12  For  this 
experiment,  the  current  necessary  to  heat  the  probe  is  provided  by  a  programmable  power 
supply  with  floating  outputs.  At  each  location  in  the  plume,  the  current  is  steadily 
mcreased  and  the  potential  with  respect  to  ground  at  the  negative  terminal  of  the  supply  is 
recorded.  This  method  allows  for  verification  of  a  well-defined  plateau  in  the  voltage- 
current  trace,  which  indicates  plasma  sheath  neutralization.  Considering  that  the  voltage 
drop  across  the  emitting  filament  never  exceeds  6  V,  the  uncertainty  in  the  plasma 
potential  measurements  is  estimated  to  be  ±3  V. 

Electrostatic  Energy  Analyzer 

A  45°,  parallel-plate  type  electrostatic  energy  analyzer  (ESA)  consists  of  two 
parallel  plates  separated  by  a  distance,  d.  One  of  the  plates  is  electrically  grounded  while 
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the  other  is  biased  to  a  positive  potential,  VR,  to  reflect  ions  admitted  through  a  slit  in  the 
grounded  baseplate.  After  being  deflected  by  the  applied  electric  field,  ions  of  a  selected 
initial  velocity  to  charge  ratio,  m/qj,  pass  through  a  second  slit  and  are  collected  by  a 
detector  as  illustrated  in  Fig.  4.  Since  only  ions  of  a  specific  energy  to  charge  ratio  are 
collected,  the  ESA  acts  as  a  velocity  per  charge  filter  for  a  given  ion  species.  For  the 
specific  case  of  a  45°  ion  injection  angle,  the  properties  of  the  collected  ions  are  related  to 
the  voltage  of  the  repelling  plate  by  Eqn  3.  Introducing  the  equivalent  ion  voltage  allows 
the  relationship  to  be  written  in  the  simple  form  of  Eqn.  4.  Thus,  the  collector  current 
measured  as  a  function  of  the  applied  plate  voltage  is  proportional  to  the  ion  energy  per 
charge  distribution.  Note  that  this  result  is  independent  of  the  ion  mass. 


Fig.  4:  Simplified  schematic  of  a  parallel  plate  electrostatic  energy  analyzer. 


mp?  f  L  \r 

(3) 

2 q,  ~\2d)* 

(4) 

The  main  body  of  the  ESA  consists  of  a  cube  constructed  of  mica  dielectric  and 
measuring  approximately  300  mm  (12”)  in  each  dimension.  It  is  very  similar  in  size  to  an 
instrument  used  successfully  by  Pollard  to  study  a  Hall  thruster  plume.13  The  parallel 
plates  are  constructed  of  1.6  mm  thick  aluminum  and  have  the  following  relevant 
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dimensions;  L-152.4  mm,  d=76.2  mm,  and  w=1.5  mm.  Two  field  correction  plates  are 
placed  between  the  main  plates  at  equal  intervals  and  biased  by  resistor  strings  to  reduce 
the  adverse  effects  of  fringing  electric  fields.  During  data  collection,  VR  is  swept  from  0 
to  600  Volts  and  the  resulting  current  to  the  grounded  collector  plate  is  recorded  using  a 
picoammeter.  Multiple  ESA  traces  recorded  at  each  data  point  demonstrate  excellent 
repeatability,  as  illustrated  elsewhere.5 
Retarding  Potential  Analyzer 

The  retarding  potential  analyzer  (RPA)  diagnostic  allows  the  collection  of 
selectively  filtered  ions  by  applying  a  retarding  potential  across  an  inlet  grid.  For  a  given 
grid  potential,  only  ions  with  energy  to  charge  ratios  greater  than  the  grid  voltage  pass 
through  the  inlet  and  reach  the  collector.  The  magnitude  of  the  derivative  of  the  resulting 
current- voltage  characteristic  is  then  proportional  to  the  ion  energy  per  charge 
distribution.  The  RPA  used  in  these  experiments  is  based  on  the  multi-gridded  energy 
analyzer  design  of  Hutchinson.14  It  is  composed  of  three  grids  and  is  shown 
schematically  in  Fig,  5. 


Grid  1  -  Floating 

-  Grid  2  -  Electron  Repelling 
-  Grid  3  -  Ion  Retarding 


Copper  Collector 
j-  Phenolic  Sleeve 

Stainless  Steel  Body 

..  \ 


pf 

] 

i 

i 

i 

1 

uu 

Macor  Insulator 
Washers 


Fig.  5:  The  retarding  potential  analyzer. 
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The  outer  body  of  the  RPA  is  constructed  of  3 1 6  stainless  steel  (SS)  tubing,  and  a 
phenolic  sleeve  placed  inside  the  body  provides  electrical  isolation  of  the  grids.  Each  of 
the  identical  grids  are  cut  from  0. 127  mm  (0.005”)  thick  3 16  SS,  photoehemically 
machined  sheet  with  0.279  mm  (0.01 1”)  diameter  openings.  The  grids  have  an  open  area 
fraction  of  38%  and  are  separated  by  Macor  washers. 

During  operation,  grid  1  is  floated  to  provide  a  non-perturbing  interface  between 
the  probe  and  the  plasma  while  grid  2  is  biased  30  V  below  ground  to  repel  electrons. 
Grid  3  is  swept  from  0  to  600  V  with  respect  to  ground  and  the  resulting  current  to  the 
collector  is  measured  using  a  picoammeter.  A  sample  RPA  current-voltage  trace,  a  cubic 
spline  fit  to  the  data,  and  the  corresponding  ion  energy  distribution  function  are  shown  in 
Fig.  6.  These  data  were  taken  0.5  m  downstream  of  the  exit  plane  of  a  single  thruster,  5° 
off  centerline. 


Fig.  6.  Sample  RPA  data  and  resulting  ion  energy  profile. 

Gaussmeter 

The  magnetic  field  downstream  of  the  cluster  is  recorded  using  an  FW  Bell  model 
7030  three-axis  gaussmeter.  All  measurements  are  recorded  without  the  thrusters  in 
operation.  Although  recent  work  has  shown  the  magnetic  field  profiles  inside  an 
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operating  Hall  thruster  to  deviate  from  the  applied  profiles  due  to  fields  induced  by  the 
azimuthal  electron  drift,15  the  difference  is  expected  to  be  negligible  for  the  low-power 
thrusters  studied  here  because  of  the  low  current  levels  involved.  The  magnetic  field 
profiles  presented  in  this  paper  are,  therefore,  believed  to  be  realistic  representations  of 
those  that  occur  downstream  of  an  operational  cluster. 

Results  and  Discussion 

Magnetic  Field 

Figure  7  shows  magnetic  field  data  recorded  in  the  XZ  plane  of  thrusters  2  and  3 
and  in  the  YZ  plane  of  thrusters  3  and  4.  The  differences  in  these  plots  are  attributable  to 
the  different  direction  of  magnet  current  flow  between  thrusters  2  and  4,  Thrusters  2  and 
3  are  operated  with  the  electromagnets  in  the  nominal  configuration  while  the  current 
flow  was  reversed  in  thruster  4.  Reversing  the  polarity  of  electromagnets  in  alternate 
thrusters  of  a  cluster  has  been  suggested  as  a  means  of  canceling  the  disturbance  torques 
that  typically  result  from  the  slight  ExB  drift  of  the  beam  ions.1,2,16  The  data  presented  in 
Figs,  6  and  7  will  be  used  to  test  the  previously  published  theory  that  the  plasma  potential 
profiles  of  a  cluster  can  be  predicted  from  magnetic  field  data.4 
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Fig.  7:  Magnetic  field  profiles  downstream  of  a.)  thrusters  2  &  3  and  b.)  thrusters  3  &4. 

Plasma  Density 

A  triple  Langmuir  probe  is  used  to  measure  the  plasma  number  density  at  5  mm 
intervals  in  the  cluster  plume.  Data  are  recorded  in  both  the  XZ  plane  of  thrusters  2  and  3 
and  the  YZ  plane  of  thrusters  3  and  4.  For  both  planes,  data  are  recorded  with  each 
thruster  operating  alone  and  with  two  thrusters  operating  simultaneously.  Due  to  the 
good  agreement  between  the  two  data  sets,  only  the  data  recorded  in  the  YZ  plane  of 
thrusters  3  and  4  are  reported  here. 

The  plasma  density  profiles  downstream  of  thrusters  3  and  4  are  shown  in  Fig.  8. 
As  this  plot  shows,  the  maximum  number  density  50  mm  downstream  of  the  cluster  exit 
plane  is  roughly  1x10  m  .  This  value  decreases  rapidly  in  the  downstream  direction 
and  by  Z— 250  mm  the  maximum  plasma  density  has  decreased  by  more  than  an  order  of 
magnitude  to  about  5x10  m3.  Figure  8  shows  a  well-defined  jet  structure  downstream 
of  each  individual  thruster.  By  about  250  mm  downstream  the  plumes  have  merged  to 
the  pomt  that  the  density  is  nearly  constant  across  the  width  of  the  cluster  and  resembles 
file  profile  that  would  be  expected  downstream  of  a  large  monolithic  thruster. 
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Fig.  8;  Plasma  density  downstream  of  thrusters  3  and  4. 

Figure  9  shows  plasma  density  profiles  at  axial  distances  of  50, 150,  and  250  mm 
downstream  of  the  cluster  exit  plane.  The  unmarked  lines  in  these  plots  are  obtained  by 
linear  superposition  of  the  data  recorded  with  thruster  3  and  thruster  4  running 
independently.  The  measurements  of  plasma  density  taken  with  both  thrusters  operating 
simultaneously  agree  with  the  calculated  values  to  well  within  the  margin  of  error  of  the 
triple  probe  diagnostic.  This  implies  that  the  density  in  a  cluster  plume  can  be  predicted 
by  summing  the  contributions  of  each  individual  thruster,  as  shown  in  Eqn.  5. 
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Fig.  9:  Plasma  density  a.)  50  mm,  b.)  150  mm,  and  c.)  250  mm  downstream  of  thrusters  3 

and  4. 

Electron  Temperature 

The  electron  temperature  contours  recorded  downstream  of  thrusters  3  and  4  are 
displayed  in  Fig.  10,  The  temperature  varies  between  roughly  3  eV  at  Z=50  mm  along 
the  thruster  centerlines  to  less  than  1  eV  near  the  boundaries  of  the  sampled  region.  The 
data  show  slight  discrepancies  in  the  electron  temperature  in  die  near-field  of  each 
individual  thruster.  Measurements  recorded  downstream  of  thrusters  2  and  3  (not  shown) 
indicate  similar  differences,  thus  the  variations  are  not  believed  to  be  a  result  of  the 
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reversed  magnetic  field  profiles  mentioned  previously.  Rather,  the  discrepancies  are 
probably  due  to  tolerances  in  the  manufacturing  process  or  differences  in  the  cumulative 
time  of  operation  between  the  devices.  The  difference  in  the  electron  temperature  in 
front  of  each  thruster  decreases  as  a  function  of  downstream  distance  and  by  roughly 
Z—90  mm  the  difference  between  the  two  units  becomes  negligible. 


Fig,  10:  Electron  temperature  profiles  downstream  of  thrusters  3  and  4. 

Electron  temperature  traces  measured  at  axial  locations  of  50,  150,  and  250  mm 
are  shown  in  Fig.  1 1.  The  unmarked  line  in  each  plot  is  calculated  using  Eqn.  6  and 
represents  a  first  attempt  to  predict  the  electron  temperature  in  the  cluster  plume.  The 
simplistic  approach  of  calculating  a  density  weighted  average,  as  indicated  by  Eqn.  6, 
seems  to  slightly  underpredict  the  measured  temperature,  particularly  in  the  region 
between  the  thrusters.  The  electron  temperatures  measured  during  thruster  operation  are 
consistently  higher  than  those  recorded  during  single  thruster  operation;  however  the 
difference  is  generally  less  than  0.2  eV,  which  is  within  the  uncertainty  of  the  diagnostic. 
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shown  in  Fig.  12.  An  interesting  feature  shown  in  this  plot  is  the  unique  plasma  potential 
profile  in  the  area  between  the  thrusters.  Between  approximately  Y=-30  and  Y=30  mm, 
the  plasma  potential  increases  with  downstream  distance  indicating  that  there  exists  a 
region  where  the  electric  field  vector  is  oriented  in  the  upstream  direction.  This  can  be 
seen  clearly  in  Fig.  13,  which  shows  the  plasma  potential  profiles  at  various  axial 
locations.  The  reversed  electric  field  could  potentially  cause  ions  produced  in  the  area 
between  the  thrusters  to  be  accelerated  upstream  toward  the  spacecraft  on  which  the 
thrusters  are  mounted.  Although  this  could  hypothetically  result  in  an  increased  erosion 
rate  in  some  areas  due  to  increased  ion  impingement,  the  effect  is  expected  to  be 
negligible  since  the  impinging  ions  are  unlikely  to  experience  accelerating  potentials 
greater  than  a  few  volts  in  the  reverse  direction. 


Fig.  12:  Plasma  potential  downstream  of  thrusters  3  and  4. 


Fig.  13:  Plasma  potential  profiles  downstream  of  two  operating  Hall  thmsters. 
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It  has  been  suggested  that  the  plasma  potential  profiles  downstream  of  a  cluster 
could  be  predicted  by  simply  integrating  the  magnetic  field  data.4  This  is  contradicted  by 
the  measurements  presented  in  Fig.  14,  which  shows  the  plasma  potential  downstream  of 
thrusters  3  and  4  at  axial  distances  of  60, 100,  and  140  mm.  Clearly,  integration  along 
the  magnetic  field  lines  depicted  in  Fig.  7b  does  not  result  in  the  observed  potential 
profiles. 


Fig.  14:  Plasma  potential  profiles  at  a.)  Z=60  mm,  b.)  Z=100  mm,  and  c.)  Z=140  mm. 

A  more  conventional  method  for  relating  the  magnetic  field  architecture  to  the 
plasma  potential  involves  consideration  of  electron  dynamics  in  a  plasma.  Along  a 
magnetic  field  line,  the  motion  of  electrons  is  governed  purely  by  electrostatic  forces  and 
can  be  described  by  the  well  known  Boltzmann  relation.17  This  leads  naturally  to  the 
definition  of  a  thermalized  potential,  fa,  which  is  given  by  Eqn.  7  and  conserved  along  a 

line  of  force.  In  the  derivation  of  Eqn.  7,  the  electron  temperature  has  been  assumed 
constant  along  lines  of  force. 
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The  concept  of  thermalized  potential  is  useful  in  the  design  of  Hall  thrusters  since 
it  shows  that  the  magnetic  field  lines  can  be  approximated  as  equipotential  lines  in 
situations  where  the  electron  temperature  is  negligible  compared  to  the  plasma  potential. 
In  other  words,  the  thermalized  potential  is  a  useful  tool  for  predicting  plasma  potential  in 
situations  where  electrons  are  tightly  bound  to  the  magnetic  field  lines.  This  method, 
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however,  is  less  useful  in  the  thruster  plume  since  the  correction  term  due  to  thermal 
effects  and  density  gradients  can  be  as  large  as  the  plasma  potential. 

Comparing  the  plasma  potential  data  of  Fig.  12  to  the  magnetic  field  profiles 
shown  in  Fig.  7b,  it  is  clear  that  the  lines  of  force  do  not  correspond  to  equipotential 
contours.  This  is  not  surprising,  since  the  magnetic  field  strength  is  generally  less  than 
10  G  and  the  electrons  are  only  weakly  magnetized  throughout  the  areas  where  the 
plasma  potential  is  presented.  In  this  situation,  thermal  effects  and  density  gradients  are 
dominant  over  the  effects  of  the  magnetic  field,  and  the  plasma  potential  is  described  by 
the  Boltzmann  relation  given  by  Eqn.  8. 19  The  profiles  calculated  using  Eqn.  8  are  shown 
in  Fig.  14  and  generally  agree  to  within  one  volt  of  the  measured  values,  except  in  the 
most  upstream  locations  of  the  sampled  region.  In  utilizing  Eqn.  8,  the  reference  density 
(n0=7xl017  m'3  in  this  case)  is  chosen  so  as  to  make  the  plasma  potential  calculated  along 
the  centerline  of  thruster  3  at  Z=  100  mm  match  the  measured  value.  While  the  choice  to 
match  the  value  at  100  mm  is  arbitrary,  a  similar  approach  is  expected  to  be  valid  in  most 
practical  cluster  configurations  since  the  data  presented  here  shows  the  plasma  potential 
directly  downstream  of  one  thruster  to  be  largely  unaffected  by  the  surrounding  devices. 
Implementation  of  Eqn.  8  along  with  Eqns.  5  and  6  thus  allows  the  most  basic  plasma 
properties  downstream  of  a  cluster  of  identical  Hall  thrusters  to  be  predicted  based  solely 
on  measurements  or  simulations  of  a  single  unit.  Results  obtained  in  this  way  appear  to 
be  accurate  to  within  the  margin  of  error  of  typical  plasma  diagnostics. 
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Single  Thruster  Ion  Energy  Profiles 


As  an  initial  test  of  the  energy  analysis  diagnostics,  energy  distributions  measured 
with  the  RPA  are  compared  to  those  measured  with  the  ESA  for  a  single  thruster.  Figure 
15  shows  the  measured  distribution  function  for  a  single  thruster  at  0°,  15°,  and  30°  off 
centerline  for  each  instrument  Notice  the  relatively  good  agreement  between  the  two 
devices.  For  each  of  the  three  angular  locations,  the  voltage  at  which  the  peak  in  the 
distribution  function  occurs  agrees  to  within  8  volts.  For  example,  on  the  thruster 

centerline  the  primary  peak  was  measured  at  220  volts  by  the  RPA  and  at  228  volts  by 
the  ESA. 


a. 


b. 
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Fig.  15:  Comparison  of  ESA  and  RPA  diagnostics  0,5  m  downstream  of  a  single  Hall 
thrusters  at  angles  of  a.)  0,  b.)  15,  and  c.)  30  degrees  off  centerline. 

The  most  noticeable  difference  between  the  diagnostics  demonstrated  by  Fig.  15 
is  the  appearance  of  secondary  peaks  at  voltages  above  and  below  the  primary  ion 
voltage,  which  are  more  pronounced  in  the  ESA  traces.  Additionally,  the  primary  peak  in 
the  distribution  is  consistently  wider  when  measured  with  the  RPA  as  opposed  to  the 
ESA.  The  shape  of  the  distribution  function  is  likely  to  be  more  accurate  in  the  ESA 
traces  since  those  data  are  not  subject  to  the  effects  of  numerical  differentiation.  The 
location  of  the  primary  peak  in  the  distribution,  however,  is  likely  to  be  more  accurately 
depicted  by  the  RPA  because  slight  misalignment  of  the  grid  components  would  not  be 
expected  to  alter  the  performance  of  this  device.  Slight  misalignment  or  improper 
spacing  of  the  plates  in  the  ESA,  on  the  other  hand,  could  cause  a  shift  of  several  volts  in 
the  measured  distributions. 

Figure  16  summarizes  the  ion  energy  distributions  recorded  by  the  ESA,  while 
Fig.  17  depicts  similar  data  recorded  by  the  RPA.  All  data  are  recorded  0.5  meters 
downstream  of  the  thruster  along  a  radial  arc.  Although  traces  have  been  recorded  for 
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both  positive  and  negative  angles  off  centerline,  only  data  for  the  positive  angles  are 
reported  here  due  to  the  high  degree  of  symmetry  exhibited  by  the  plume.  The  ESA 
traces  show  the  peak  ion  energy  to  charge  ratio  to  occur  at  approximately  228  volts  for 
most  of  the  angular  spectrum,  while  the  RPA  shows  the  peak  at  220  volts.  The  secondary 
structure  occurring  at  energy  to  charge  ratios  below  150  volts  can  be  attributed  to 
elastically  scattered  primary  ions.13, 19  The  high-energy  population  shown  at  voltages  in 
excess  of  the  discharge  voltage,  particularly  at  low  angles  off  centerline,  is  likely  due  to 
beam  ions  that  have  undergone  charge  decreasing  collisions.13,20  Data  are  not  shown  for 
the  ESA  at  angles  greater  than  60°  due  to  the  prohibitively  small  signal  to  noise  ratio  in 
this  regime.  RPA  data,  however,  show  the  plume  to  be  composed  primarily  of  low- 
energy  charge  exchange  products  at  angles  greater  than  70°. 
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Fig,  17:  Single  thruster  RPA  data  at  angles  of  a.)  0-25,  b.)  30-60,  c.)  70-90  degrees  off 

centerline. 


Cluster  Ion  Energy  Profiles 

Unlike  the  data  recorded  for  the  single  thruster  case,  the  measurements  obtained 
with  the  ESA  and  RPA  aligned  to  the  center  of  the  cluster  show  marked  differences 
between  the  two  diagnostics.  These  differences  are  believed  to  be  caused  primarily  by 
the  different  acceptance  angles  of  the  RPA  and  ESA.  The  ESA  entrance  slit  provides  an 
ion  acceptance  angle  of  approximately  4°  in  one  direction  and  0.5°  in  the  other  direction, 
while  the  cylindrical  RPA  has  an  acceptance  cone  half  angle  of  approximately  25°.  This 
discrepancy  is  not  important  for  the  case  of  a  single  thruster  because  both  diagnostics  are 
able  to  image  the  entire  width  of  the  thruster  at  a  downstream  distance  of  0.5  m.  At  this 
distance,  the  ESA  images  a  cross  section  only  about  70  mm  wide.  In  the  cluster 
configuration,  this  results  in  the  ESA  imaging  the  space  between  the  thrusters  rather  than 
the  thrusters  themselves.  The  RPA,  on  the  other  hand,  has  a  sufficiently  wide  viewing 
angle  to  accept  ions  originating  from  any  of  the  four  thrusters. 
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Figure  18  summarizes  the  cluster  data  collected  with  the  parallel  plate  energy 
analyzer.  At  angular  positions  less  than  10°  with  respect  to  centerline,  the  peak  in  the 
distribution  occurs  at  energy  to  charge  ratios  near  the  250  volt  discharge  voltage. 
Between  10°  and  20°  the  peak  shifts  down  to  approximately  134  volts,  which  is  near  the 
voltage  of  the  elastically  scattered  ions  measured  in  the  plume  of  a  single  thruster.  The 
134  volt  peak  can  be  observed  out  to  80°  off  the  cluster  axis  before  the  signal  is  lost 
between  80°  and  90°.  It  should  be  noted  that  the  signal  level  recorded  by  the  ESA  in  this 
configuration  is  approximately  a  factor  of  25  lower  than  that  measured  for  a  single 
thruster  due  to  the  limited  viewing  angle. 


a. 


b. 
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Fig.  18:  Ion  energy  profiles  downstream  of  the  cluster  at  angles  of  a.)  0-15  and  b.)  15-80 

degrees  off  centerline. 


The  RPA  data  presented  in  Fig.  19  show  several  unusual  characteristics, 
particularly  along  the  cluster  centerline  where  the  spectrum  shows  three  distinct, 
repeatable  peaks  at  224, 116,  and  74  volts.  As  explained  by  Gallimore,  the  peaks  at  1 16 
and  74  volts  could  be  caused  by  ions  exiting  the  thruster  at  a  beam  kinetic  energy,  Vb,  of 
224  volts  before  undergoing  charge  exchange  (CEX)  collisions  that  result  in  populations 
with  energy  to  charge  ratios  of  approximately  Vb/2  and  Vj/3,  respectively.20  Just  5°  off 
centerline,  however,  the  spectrum  changes  to  a  double  peaked  structure  with  equally 
abundant  populations  occurring  at  122  and  222  volts.  As  the  angle  off  centerline  is 
increased,  the  two  peaks  merge  together  to  form  a  single  peak  near  206  volts  with  a  low 
energy  tail  as  shown  in  Fig.  19b.  Charge  exchange  collisions  tend  to  generate  signatures 
in  the  energy  per  charge  spectra  at  discrete  multiples  of  the  main  peak  in  the 
distribution.  The  fact  that  the  ratio  between  the  observed  peak  voltages  changes  as  a 
function  of  angle  seems  to  indicate  that,  if  the  multi-peak  structure  is  due  to  collisions, 
the  dominant  interactions  are  likely  elastic  rather  than  CEX.  The  existence  of  elastic 
scattering  peaks  at  voltages  below  the  main  peak  is  consistent  with  numerical  simulations 
that  take  into  account  recently  calculated  cross  sections  for  collisions  between  xenon  ions 
and  neutrals.  The  shift  of  the  low  energy  population  to  higher  voltage  with  increasing 

angle  seen  in  Fig.  19a,  however,  is  opposite  to  previously  observed  trends.13  This 
suggests  that  a  physical  mechanism  other  than  collisions  may  contribute  to  the  low- 
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energy  structure.  At  angles  greater  than  50°  with  respect  to  the  thrust  axis,  the  spectra  are 
dominated  by  a  low-energy  population  that  shifts  to  lower  voltage  with  increasing  angles. 


a. 


b. 
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Fig.  19:  RPA  data  at  a.)  low,  b.)  medium,  and  c.)  high  angles  off  the  cluster  centerline. 

In  addition  to  collisions,  a  possible  contributing  factor  to  the  low-energy  structure 
involves  ion  focusing  as  a  result  of  clustering.  As  shown  in  Figs.  12-14,  the  plasma 
potential  profile  downstream  of  a  cluster  is  fundamentally  different  than  the  profile  in  the 
plume  of  a  single  thruster.  When  ions  exit  a  single  Hall  thruster,  they  experience  a 
continuous  decline  in  plasma  potential  regardless  of  the  direction  in  which  they  exit  the 
thruster.  In  other  words,  the  electric  field  vector  is  everywhere  directed  away  from  the 
device.  When  multiple  thrusters  are  operated  together,  however,  a  minimum  in  the 
plasma  potential  occurs  in  the  region  between  the  thrusters.  This  results  in  a  situation 
where  an  ion  directed  toward  the  center  of  the  cluster  can  be  deflected  downstream  by  the 
plasma  potential  hill  created  by  adjacent  thrusters.  This  situation  is  sketched  in  Fig.  20 
below,  in  which  the  dashed  lines  represent  contours  of  constant  plasma  potential  and  the 
heavy  arrows  represent  the  paths  of  sample  ions.  The  phenomenon  illustrated  in  Fig.  20 
may  lead  to  ion  focusing  in  which  ions  initially  directed  toward  the  cluster  center  are 
deflected  to  lower  angles  with  respect  to  the  cluster  centerline.  This  effect  may  be 
responsible  for  the  slightly  reduced  beam  divergence  reported  by  Hargus  and  Reed  for 
two  operating  thrusters  compared  to  that  predicted  by  linear  superposition  of  the  ion  flux 
from  individual  thrusters.4 


Fig,  20:  Ion  focusing  mechanism  due  to  clustering. 

The  mechanism  of  ion  focusing  presented  here  requires  additional  analysis  before 
its  effects  on  the  plasma  plume  can  be  evaluated  quantitatively.  However,  one  can  gain 
insight  into  several  aspects  of  the  energy  spectrum  by  resorting  to  a  simple 
phenomenological  discussion.  Consider  two  ions,  A  and  B,  exiting  a  thruster  and 
traveling  in  an  identical  direction  toward  the  center  of  the  cluster,  but  with  different 
initial  kinetic  energies.  In  this  situation,  the  slower  moving  ion,  B,  would  be  deflected  by 
a  given  potential  rise  to  a  greater  extent  than  its  high-energy  counterpart,  ion  A,  as 
depicted  in  Fig.  20.  Considering  this,  a  detector  swept  through  the  plume  would  detect 
ion  A  at  a  higher  angle  off  centerline,  while  ion  B  with  its  lower  energy  would  be 
deflected  further  downstream  and  detected  at  a  relatively  low  angle.  This  phenomenon 
may  account  for  the  secondary  structure  shown  in  Fig,  19a,  in  which  the  low  energy 
population  shifts  to  higher  voltages  with  increasing  angle  off  centerline.  This  feature  is 
only  observable  during  cluster  operation  and  is  most  pronounced  at  low  angles. 
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In  addition  to  the  cluster  configuration  discussed  above,  data  have  also  been 
obtained  with  four  thrusters  in  operation,  but  with  the  ESA  and  RPA  aligned  to  the 
centerline  of  TH  3  rather  than  die  centerline  of  the  cluster.  These  data  are  presented  in 
detail  elsewhere5  and  are  qualitatively  similar  to  the  results  discussed  above.  The  most 
notable  features  observed  in  this  “offset  cluster”  configuration  are  die  tendency  of  the 
primary  peak  in  the  distribution  to  shift  to  lower  voltages  with  increasing  angle  off 
centerline  and  the  predominance  of  structures  at  energies  well  below  the  discharge 
voltage  at  angles  greater  than  40°,  Although  the  secondary  structure  is  visible  in  the 
single  thruster  data,  its  magnitude  is  much  larger  during  cluster  operation,  which  seems  to 
indicate  a  dramatic  increase  in  ion  scattering  due  to  multi-thruster  operation. 

Conclusions 

A  combination  of  triple  Langmuir  probes  and  floating  emissive  probes  is  used  to 
characterize  the  electron  number  density,  electron  temperature,  and  plasma  potential  in 
the  plume  of  a  low-power  Hall  thruster  cluster.  The  results  show  that  properties  in  the 
cluster  plume  can  be  predicted  using  simple  analytical  relations  and  knowledge  of  a 
single  thruster  plume.  In  particular,  the  plasma  density  can  be  predicted  by  simple  linear 
superposition  while  the  electron  temperature  is  estimated  to  be  a  weighted  average  of  that 
due  to  the  plumes  of  individual  thrusters.  Having  obtained  predictions  of  electron 
number  density  and  electron  temperature,  these  results  are  then  used  to  estimate  the 
resulting  plasma  potential  via  the  Boltzmann  relation. 

Ion  energy  spectra  are  obtained  using  both  a  parallel  plate  electrostatic  analyzer 
and  a  retarding  potential  analyzer  at  locations  0.5  meters  downstream  of  the  thruster 
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cluster.  Compared  to  results  obtained  in  the  plume  of  a  single  thruster,  these  data 
indicate  a  profound  increase  in  the  fraction  of  low-energy  ions  measured  at  voltages 
below  the  primary  peak  in  the  distribution.  This  feature  is  believed  to  be  due  to  an 
increase  in  elastic  scattering  of  beam  ions  as  well  as  the  unique  plasma  potential  profiles 
downstream  of  a  cluster.  Ion  focusing  as  a  result  of  the  plasma  potential  distribution  is 
hypothesized  to  be  responsible  for  the  response  of  low-energy  structures  reported  here,  as 
well  as  effects  observed  in  recently  published  ion  flux  data. 
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